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Abstract 

In this work we apply Wang-Landau simulations to a simple model which has exact solutions 
both in the microcanonical and canonical formalisms. Lhe simulations were carried out by using an 
updated version of the Wang-Landau sampling. We consider a homopolymer chain consisting of N 
monomers units which may assume any configuration on the two-dimensional lattice. By imposing 
constraints to the moves of the polymers we obtain three different models. Our results show that 
updating the density of states only after every N monomers moves leads to a better precision. We 
obtain the specific heat and the end-to-end distance per monomer and test the precision of our 
simulations comparing the location of the maximum of the specific heat with the exact results for 
the three types of walks. 
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I. INTRODUCTION 



The sequencing of the Human Reference Genome, announced ten years ago, provided a 
roadmap that is the basic foundation for modern biomedical research 1] . This monumental 
achievement was enabled by developments in DNA (homopolymer) sequencing technology 
that allowed data production which exceed the original descripion of Sanger sequencing |2j. 
Linear polymers are the simplest physical systems that can be studied in the framework of 
random walks models. They are long chain-like molecules formed by repetition of a basic 
unit or segment, where more importantly the polymer is flexible, i.e., it can assume different 
geometric configurations. 

Recently, the study of homopolymers has been established by various techniques in con- 
densed matter physics. Cohen et al. js], studied the behavior of single file translocation of 
a homopolymer through an active channel under the presence of a driving force by using 
Langevin dynamics simulation. Previous works on homopolymers which studied the denat- 
uration of circular DNA are extensions of the Poland-Scheraga model ^]. Experimentally, 
the viscoelastic properties of a binary mixture of a mesogenic side-chain block copolymer in 



or mass concentrations ranging 

a. 



a low molecular weight nematic liquid crystal are studied 
from the diluted regime up to a liquid crystalline gel state 

Although Monte Carlo simulations play an important role for the study of phase tran- 
sitions and critical phenomena, some well-known difficulties arise when one uses standard 

fl 

algorithms (one-flip algorithms) [6] for the study of random walks models. These difficul- 
ties have been overcomed by the development of alternative Monte Carlo methods, such as 



parallel-tempering [7], cluster algorithms [8|, multicanonical algorithms [9|, and more re 



cently the Wang-Landau method 



[9] 



10] . This method has been applied with great success to 



11 



13|- 



many systems, in particular to polymers in lattice 

In the present paper, using the Wang-Landau method, we investigate the computer sim- 
ulations of a homopolymer model with exact solution in canonical and microcanonical for- 
malism. In section [III we give an introduction to the three studied models and we briefly 
present the mathematical background. In section UTTl we outline shortly how the simulations 
were carried out. In section [TV] we show and discuss the results for all the models. 
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II. MODELS AND FORMALISM 



A. Model 1 

We consider a homopolymer chain of N monomers units of length a which may assume 
any configuration on a two-dimensional latice [l4j . One end of the polymer is fixed and its 
taken as the origin of coordinates, shown by an open circle in Figure HJ a). The other end 
of this linear chain is subject to an externally applied tension r, acting along the positive 
x-axis. A possible realization of this model is sketched in Figure HJ a). Each polymer unit 
is permitted to lie either parallel or antiparallel to the x-axis and we assign the works —ra 
and +ra to these two orientations. We denote by L x the distance between the ends of the 
polymer chain, on the x-direction. Each monomer unit has the additional possibility of lying 
perpendicular to the x-axis, in the +y or — y directions. We associate a positive energy e 
to such a perpendicular monomer and the distance between the ends of the chain on the 
y-direction is denoted by L y . 

The hamiltonian of this model can be written as 



where N£ and N x are the number of monomers along the +x and — x directions respectively, 
and similarly for iV+ and N~. Since the tension in the y direction is zero, we can assume 



Hi = {N+ + N;) e + (N- - N+) ra 



(1) 



N+ = N~. Then 



n = n: + n- + n+ + n, 



y ' 



(2) 



K-N- = L 



(3) 



N+ + N, 



V 



U. 



(4) 



from which we find 




(5) 



N~ = -(N-U-L X ), 



(6) 
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and 



K = N v = \U. (7) 

The number of configurations of the polymer consistent with a given end-to-end distance 
in the x-direction, L x (the dimensionless length of the polymer), and U (the number of 
monomers lying perpendicular to the x-axis), is 

An 

n (L x , U, AO = ■ . (8) 

From Eq. (0Q), if we set e = ra = 1, we can write for a given energy level 

E = U-N+ + N~, (9) 

and from Eq. fl2]) 

N = U + N+ + N~. (10) 
Adding the equations (jSJ) and ( TTUjl we obtain 

U + K = ^. (ID 



Inserting Eqs. (Ej), (EJ) and (CD into (JH1) , using Eqs. (Ej), (JTUJ), and ([TT)), and setting 
A^~ = n we obtain the number of configurations with energy E as 

JV+E 

m = hh¥¥W¥^M^-' (12) 

Using the definition of the entropy and defining dS/dU = 1/T, we obtain [14]: 

L x N sinh (ra/ksT) 

N cosh (ra/ksT)! + exp (— e/k B T) 

B. Model 2 



(13) 



In this model we consider that each polymer unit is parallel to the x-axis; no antiparallel 
move is allowed. Additionally the polymer units have the possibility of lying on the +y and 
— y directions. In Figure [T] b) is sketched a possible configuration of this model, where we 
depicted by an impenetrable wall the forbidden region along the 
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The hamiltonian for the model 2 can be written as 



U 2 = (N+ + N y ) e — N+tcl, (14) 

where N x is the number of monomers along the +x direction and similarly for Ny and N~ . 
Assuming again Ny = N~ we obtain the following equations: 

NZ= l -{N-U + L x ), (15) 

and 

K = K = % (is) 

For this model the number of configurations of the polymer consistent with a given L x 
and U is 

n (L x , U, N) = : . (17) 

x y y ' 

From Eq. f|T4l . if we set £ = ra = 1, we can write the energy as 

E = U-N+, (18) 
and for N = N+ + N~ + iV+ we obtain 

N = U + N+. (19) 

After similar calculations as for model 1, subsection III Al we obtain the density state and 
end-to-end distance per monomer, respectively 

N+E 

TV' 



n=0 

and 

U 1 



V l + 2ex P (-g|) 



(21) 
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C. Model 3 



In this model we consider that each polymer unit is parallel to the x-axis and we allow 
only the possibility of lying on the +y direction. Thus, we restrict to the situation of positive 
values for both axes. In Figure [T]c) we show a possible configuration of this model. 

The hamiltonian for the model can be written as 

U 3 = N+e - N+ra, (22) 

where A+ is the number of monomers along the +x direction and similarly for iV+ . Then 
Using the method previously described, we obtain N£ and for this model 

Nt = \{N-U + L x ) (23) 



and 



N;= l -{N + U-L x ). (24) 



For this model the number of configurations with energy E is given by 

N 

Using the equation for g(E) we obtain the length of the polymer as 

L 



J x 



^ _ l + exp(-g|)' <26) 



III. SIMULATIONS 



We define a Monte Carlo 



In our simulations we followed the prescriptions of Ref. 
step (MCS) as giving sequentially to any unit the possibility of changing its direction with 
identical probability to any allowed direction or remaining in the same one. At the beginning 
of the simulation we set S(E) = for all energy levels, where S(E) = \ng(E). The random 
walk in the energy space runs through all energy levels from E min to E max with a probability 

p{E^E') = minjexp (s(E)-S(E'^j (27) 

where E and E' are the energies of the current and the new possible configurations. After 
N trial moves we update H(E) -»■ H(E) + 1 and S(E) -»■ S(E) + F h where F { = In /<, 



fo = e = 2.71828... and fi + ± = yffi (where fi is the so-called modification factor and 
H(E) is a histogram accumulated for each fi). The flatness of the histogram is checked 
after a number of Monte Carlo (MC) steps and usually the histogram is considered flat if 
H(E) > 0.8(H), for all energies, where (H) is an average over the energies. If the flatness 
condition is fulfilled we update the modification factor to a finer one and reset the histogram 
H(E) = 0. The simulations are continued up to f final = /u and the microcanonical averages 
were accumulated from the very beginning (f micro = fo), results obtained by Ferrera and 



Caparica 161 ] . Having in hand the density of states, one can calculate the canonical average 



of any thermodynamic variable as 

{X)T - E E 9(E)e-^ ' (28) 
where (X)e is the microcanonical average accumulated during the simulations and /3 = 
l/kgT, is the Boltzmann constant and T is the temperature. 



IV. RESULTS AND DISCUSSION 

In Figure[T]we depicted the three models which were presented in section[TTl The difference 
between the models is given by the allowed moves. In the first model, denoted by a) in the 
figure, we allow the polymer unit to move along the positive or the negative directions of 
the x-axis, the same for the y-axis. The end-to-end distance of the polymer on x-direction 
is denoted by L x and with L y we denote this distance on y-direction. In the second model, 
b) in figure [I], the antiparallel motion along the x-axis is forbidden. In the third considered 
model there are allowed moves only on the positive side of both x and y-axis. 

Using the simulated and the exact density of states in Eq. (f2"8"j) we calculate the specific 
heat given by 

and the mean end-to-end distance 

(L x ) = (\xn — Xi\), (30) 

where E is the energy of the configurations and x\ and xn are the corresponding x- 
coordinates of the ends of the polymer. 
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In Figure [2] we plot in semi-logarithmical scale the density of states g(E) for all the three 
models. Here we plot the results obtained from the simulations (symbols in the figure) and 
also the exact theoretical results (continuous lines in the figure) for a polymer of iV = 500 
monomers. In Figure E] we plot the end-to-end distance along the x-axis, L x , as a function 
of the temperature. Here we rescale L x by the total number of monomers, N = 500 in this 
case. We observed a very good agreement between the simulation results (symbols in the 
figure) and the theoretical results (continuous lines in the figure), given by equations ([TBI . 
(l2~Tj) . and (|26|) . corresponding to model 1, model 2, and model 3 respectively. In Figure 0] 
we plot the specific heat per monomers for polymers with N = 500. The specific heat have 
a tail proportional to 1/T 2 in the high temperature limit. 

Finally, in Table I we present the location of peak of the specific heat for each model 
obtained by the Wang-Landau simulations and compare with the results calculated with the 
exact density of states. The simulations were carried out for 100 independent runs, adopting 
the 80% flatness criterion. One can see that in the three cases the exact results fall into the 
error bars. 



Case 


exact 


our results 


Model 1 


0.70299027 


0.7043(23) 


Model 2 


0.75335362 


0.7513(23) 


Model 3 


0.83180562 


0.8314(23) 



Table I: Temperatures of the peak of the specific heat from simulations, compared with 

the exact values. 



V. CONCLUSIONS 



We have carried out Wang-Landau simulations of a simple polymer model which has 
exact solutions in both the microcanonical and the canonical ensembles. Here we considered 
three two-dimensional models: in the first model we allowed moves in all possible directions, 
in the second model the moves along the negative x-axis are forbidden, while in the last 
model we allowed only moves along the positive direction for both x and y- axis. We have 
shown that updating the density of states only after each N trial moves and halting the 



simulations when fji na i = fu [16|, defined during the simulations, we obtain quite accurate 
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results, compared with the available analytical exact results. We have obtained a very- 
good agreement between the simulations and the exact results also for the studied physical 
quantities: the end-to-end distance and the specific heat. As expected, due to the difference 
of the allowed directions of motion, in the limit of high temperatures the end-to-end distance 
has the highest value for the third model and the lowest for the first model. 
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FIG. 1: Here we consider three types of walks: (a) Model 1: unrestricted motion along the x-axis 
and y-axis, (b) Model 2: is forbidden the backward motion on the x-axis, and (c) Model 3: the 
random walk moves only along the positive directions of x-axis and y-axis. 
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FIG. 2: Exact density of states g(E) for polymer size N = 500 for (a) Model 1; (b) Model 2 and 
(c) Model 3. The lines are the exacts results and the symbols are simulational results. 
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FIG. 3: End-to-end distance per monomer for polymer size N = 500 for (a) Model 1; (b) Model 2 
and (c) Model 3. The line is the exact result and the dots are the simulational results using the 
procedure for Wang-Landau simulation The error bars are less than symbols. 
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FIG. 4: Specific heat per monomer for polymer size N = 500 for (a) Model 1; (b) Model 2 and (c) 

Model 3. The line is the exact result and the dots are the simulational results using the procedure 

I 

for Wang-Landau simulation [lg]. The error bars are less than symbols. 
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